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THE POTASSIUM-ARGON METHOD GEOLOGICAL AGE 
DETERMINATION! 


ABSTRACT 


The development the method for dating minerals 
discussed. The decay scheme for radioactive potassium-40 examined. The results 
recent counting experiments have been averaged and give mean value 
29.4 2.7 beta emissions per second per gram potassium. Averaging the 
previous determinations the rate gamma decay gives branching ratio 
0.090 0.038 with the above beta decay rate. This close agreement with 
branching ratio 0.089 suggested previous paper the basis two 
determinations. 

Methods measuring the radiogenic argon content potassium minerals 
are discussed and the method use this laboratory described detail. 
Yield runs and grain size experiments have been carried out and are reported. 
ages have been calculated for four perthites and two lepidolites. 
These ages are excellent agreement with the best available ages obtained 
other methods. 


INTRODUCTION 


Two years after the discovery Thomson that alloy 
spontaneously emitted negative particles, Campbell and Wood conclusively 
proved the radioactivity potassium (46). Particularly since the proof 
its decay argon Thompson and Rowlands (47), geophysicists have been 
keenly aware the great possibilities method age 
determination. Such method would permit dating wide variety 
minerals and rocks which contain appreciable amounts potassium, and 
would permit more detailed studies geological history than have ever been 
possible before. 

Both calcium-40 and argon-40 are produced the decay potassium-40. 
But calcium very common constituent rocks and minerals, and calcium- 
forms 97% common calcium. result very difficult distinguish 
small amount radiogenic calcium from the larger amount common 
calcium usually present, and age-scale would 
limited use. However such objection applies the 
method because the amount common argon contained minerals very 
small. Moreover argon gas and its measurement small quantities 
gm.) presents fewer difficulties than would encountered 
analyzing for small quantities calcium. 
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Unfortunately the establishment argon method dating has been 
complicated number practical problems. The quantitative extraction 
and purification the radiogenic argon extremely difficult. The first attempt 
such measurement was made Harteck and Suess (19) who failed find 
argon dated potassium minerals. Farrar and Cady (9) reported the results 
number determinations which showed relationship between the 
amounts potassium and argon and the age minerals. The existence such 
relationship was first demonstrated Aldrich and Nier (2), who used their 
results find approximate values for the decay constants potassium-40. 
Mousuf first showed the very high radiogenic purity argon found 
old potassium feldspars. 

Another difficulty has been the fact that the decay constants potassium 
have not been precisely determined physical measurements. Until this 
done will impossible prove conclusively that small quantities argon 
are not lost diffusion from the minerals and potassium—argon time-scale 
cannot put independent basis. 

Despite these and other difficulties the future the 
method has never looked brighter. The great differences that once existed 
between geophysical and physical values for the potassium-40 decay constants 
are rapidly disappearing, and provisional ages have been 
reasonable accord with those obtained other methods. Within the past 
four years radiogenic argon measurements have been reported Smits and 
Gentner (42), Gerling and Pavlova (14), Fritze and Strassmann (12), Gentner, 
Prag, and Smits (13), Russell, Shillibeer, Farquhar, and Mousuf (37), Wasser- 
burg and Hayden (48), and Shillibeer, Russell, Farquhar, and Jones (40). 

The purpose this paper review the present status potassium—argon 
dating and report the results new measurements which strongly support 
this new method determining geological ages. 

THE RADIOACTIVITY POTASSIUM 


Geological ages based any radioactive decay demand accurate 
knowledge the decay rates the element concerned. Potassium-40 decays 
both beta emission and orbital electron capture calcium and argon 
respectively. knowledge two decay constants necessary predict the 
rate accumulation either calcium argon. convenient use the 
branching ratio which the ratio orbital electron captures beta 
emissions, and the decay constant for beta emission Using these decay 
constants, the accumulation radiogenic argon potassium mineral 
given 

directly proportional the number beta particles emitted from one 
gram potassium unit time. This has been measured many times 
number workers. The more recent determinations, which are listed 
Table are reasonable agreement. The mean these results gives beta 
decay constant 0.503 with standard deviation approxi- 
mately nine per cent. 


i 
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TABLE 
RATE BETA EMISSION POTASSIUM 


Year Authority gm. 

1948 Borst and Floyd (4) 
Graf (16) 
Hirzel and (21) 

1949 Floyd and Borst (11) 
Stout (44) 30.6 2.0 

1950 Faust (10) 
Houtermans, Haxel, and Heintz (24) 27.1 
Sawyer and Wiedenbeck (39) (ordinary potassium) 
Sawyer and Wiedenbeck (39) (enriched 0.4% 
Smaller, May, and Freedman (41) 23* 
Spiers (43) 30.5 

1951 Delaney (8) 32.0 
Good (15) 27.4 


Mean standard deviation 29.4 2.7 


*These results were omitted calculating the mean. 


The branching ratio can determined dividing the rate orbital 
capture this rate beta emission. theory the rate capture can 
most directly determined counting the number X-rays and Auger 
electrons, one which emitted the argon atom formed during 
capture process. practice both X-rays and Auger electrons are extremely 
difficult count, and precise measurements have been made. 

However there seems good reason believe that the rate orbital 
electron capture equal the rate gamma-ray emission. This supported 
precision mass measurements Johnson (26), which show that only one 
gamma ray can associated with electron capture and that none can 
associated with beta emission. has also been supported theoretical 
grounds Morrison (30). These recent results suggest that the branching 
ratio can estimated directly from the results beta particle and gamma ray 
counting experiments, suggested Bramley (5) nearly twenty years 
ago. 

Although these gamma rays are much more easily counted than argon 
X-rays, the results from gamma ray counting experiments have not been 
very good agreement. The values for the branching ratios estimated this 
way are summarized Table II. All have been calculated with the average 
rate beta decay obtained Table Where the same author authors 
have reported more than one result, the more recent value has been used; 
otherwise the list complete the best our knowledge. The average 
value for the branching ratio calculated from these measurements 0.090 
with standard deviation +0.038. obvious that results could 
selected give almost any branching ratio within this range. 

Two branching ratio measurements deserve special comment. The first 
the Auger electron counting experiment Sawyer and Wiedenbeck (39). 
This experiment more than any other indicated the equality between the 
branching ratio and the ratio gamma rays beta particles. However the 
experimental error assigned the original authors (+30%) large enough 
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TABLE 
RATE GAMMA RAY EMISSION POTASSIUM 


Quantity measured Branching 
Year Authority ratio 
radiation 
(curies/gm.K) gm. 
1930 (28) 0.061 
Miihloff (32) 3.3 0.040 
1948 Ahrens and Evans (1) 3.3 0.11 (2) 
Graf (17) 3.4+0.8 0.11(5) 
Hess and Roll (20) 8.9 0.10 (9) 
and Waffler 0.087 .012 0.087 .012 
1949 Floyd and Borst (11) 0.05 .01 0.05 
Sawyer and Wieden- 
beck (38) 2.9+0.3* 0.098 .010 
1950 Faust (10) 0.122 .014 
Houtermans, Haxel, 
and Heintz (24) 3.1+0.3 0.105 .010 
Smaller, May, and 
Freedman (41) 0.05 0.05 .01 
Spiers (43) 3.0 0.10 (2) 
Burch (6) 0.119 .006 


Mean standard deviation 0.090 0.038 


*This value has been recalculated from the original data using more recent value for the gamma 
activity potassium-42 (27). 


cause considerable uncertainty ages. The second deter- 
mination the radiogenic calcium and argon measurements made sample 
Stassfurt sylvite Inghram, Brown, Patterson, and Hess (25). They 
carried out the experiment measuring both the radiogenic calcium and 
argon with mass spectrometer isotopic dilution and obtained value 
0.126 .002 for the branching ratio. Although this experiment seems very 
good there appear two inconsistencies. First, the figures obtained 
Inghram’s group indicate age million years for the sylvite, 
approximately 100% potassium chloride. However the Stassfurt deposits are 
classed geologically mid- late-Permian (29). The age the Permian period 
has been well established 200 220 million years (22), more than twice 
the age obtained Inghram. Since this experiment the branching ratio 
determined the argon measurement, the discrepancy the age cannot 
removed assuming that the radiogenic argon was measured incorrectly, 
but can explained assuming that the mineral was not Permian age, 
that the radiogenic calcium figure was too low. the radiogenic calcium 
figure was too low the branching ratio obtained would too large. 

The second inconsistency shown studies Gentner, Prag, and Smits 
(13) which indicate that argon diffuses from sylvites with diffusion constant 
sq. cm. per sec. This value was calculated using branching ratio 
0.13. Unfortunately there record the grain size the sylvite used 
Inghram, but grain size mm. assumed, the mean size Gentner’s 
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samples, loss 42% the argon would expected. the grain size 
taken mm., the maximum size reported Gentner, the loss would 
28%. much argon had been lost from Inghram’s sample the branching 
ratio calculated from his figures would impossibly large. 

sufficient suitable samples can found, this method may eventually 
prove best, but the result obtained Inghram, Brown, Patterson, 
and Hess will some doubt until these inconsistencies are removed. 


METHODS POTASSIUM-ARGON DATING 
For calculating ages convenient use the decay 


equation the form 
40 


represents the mass radiogenic argon-40 present per gram mineral 
sample. general the gases present the mineral, including argon, are 
extracted from fusing high vacuum system. The major part the 
impurities are then removed. the purification carried out that only 
argon remains, then the volume argon may measured means 
McLeod gauge. Alternatively, known volume argon-38 (or argon-36) 
may added the unpurified gases from the sample and the volume the 
radiogenic argon calculated from mass spectrometer measurements. any 
case, mass spectrometer analysis necessary order correct for the 
presence non-radiogenic argon. The present authors use the McLeod gauge 
method, the details which will described below. 

represents the mass potassium-40 per gram sample. Nier (34) 
found that this isotope forms 0.0119% common potassium, which value 
has been substantiated (35) who obtained 0.0118%. Nier’s 
result has been used for all our calculations. The total potassium content 
determined chemical analysis. 

the partial decay constant for beta emission, has been discussed earlier 
this paper. The value 0.503 .046 obtained from the 
measurements Table has been used for all calculations. 

The value for the branching ratio, obtained averaging the results 
all gamma ray counting experiments (Table II) was 0.090 .038. value 
0.088 was obtained combining measurements reported 
Wasserburg and Hayden (48) for Besner Mine, Ontario, with lead-ratio 
age determination carried out Toronto uraninite found with the 
feldspar. The latter paper also reported calculation for dated feldspar from 
Lee Lake, Saskatchewan, which gave branching ratio 0.090. view 
these results have chosen tentative branching ratio 0.089 for our cal- 
culations. 


EXPERIMENTAL PROCEDURE 


Potassium Analyses 
Some the potassium analyses used this paper were carried out 
flame photometer Mr. Inman the Department Mineral Dressing 
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and Process Metallurgy, Ottawa. Mr. Inman has checked the accuracy 
his method using standard feldspars known potassium content obtained 
from the U.S. National Bureau Standards. Other analyses have been made 
flame photometer recently constructed Toronto Sturm. 
Several comparisons have been made Sturm the same samples analyzed 
Inman. Some these are given Table III. present consider that 


TABLE III 
INTERCOMPARISON POTASSIUM ANALYSES 


Sample No. 
Ottawa Toronto 
1069 13.07 13.05 
1063 10.96 
1025 4.38 4.07 
1059 12.25 12.01 


our potassium analyses are error not more than 7%. Since any errors 
these analyses result errors the ages, are endeavoring decrease this 
uncertainty. 

The samples from Varala, Finland, were analyzed chemically Mr. 
Wiik the Geological Survey Finland. 


The Extraction and Purification the Radiogenic Argon 


sample the mineral was finely crushed and placed container with 
least double its weight sodium hydroxide flux. The container was then 
bolted vacuum line and evacuated pressure approximately mm. 
mercury. The sample container was then disconnected from the pumping 
system means stopcock, and was connected through cold trap 
U-trap containing degassed activated charcoal. Both the cold trap and the 
U-trap were kept the temperature liquid air. liquid air temperature 
and low pressure activated charcoal adsorbs more than 99.9% the small 
volume argon present. argon condenses this temperature and pressure 
without the presence activated charcoal. This property activated charcoal 
enables use extraction line which the argon can transferred 
any desired section simply surrounding charcoal trap with Dewar 
flask containing liquid air. The main purpose the cold trap was remove 
the water released the reaction sodium hydroxide with the mineral and 
prevent the water from reaching the charcoal trap. 

The sample was heated approximately 850°C. which temperature 
the mixture was completely liquid and could stirred with metal rod. 
was maintained this temperature for hr. Water and some other im- 
purities were removed the cold trap while all the argon was adsorbed 
the activated charcoal the U-trap section. The gases the U-trap section 
were then purified with calcium furnace. The calcium was heated means 
induction heater about 800°C., which temperature calcium vapor 
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reacted rapidly with all gases, except noble gases, form solid products 
which were deposited the walls the calcium furnace. 

After purification the volume the remaining gas was measured 
McLeod gauge. ensure that purification was complete the sample was 
repurified and then remeasured. ensure that none the other noble gases 
were present photograph the line spectrum one the samples was 
taken. This showed only the line spectrum argon and mercury, the mercury 
lines being due the free mercury used the McLeod gauge. 

The sample was then adsorbed charcoal break seal tube and the tube 
cut off that the radiogenic purity the argon could measured mass 
spectrometrically. measure the relative abundances argon-40 and 
argon-36 atoms contained the sample enables one make correction 
for any contamination the sample atmospheric argon. This correction 
was most cases less than 10% and often less than 5%. Most this contami- 
nation probably due insufficient pumping the sample the beginning 
run but part may argon trapped the mineral the time 
formation. Several samples were also checked the mass spectrometer for 
hydrogen, helium, and krypton and none these gases was ever found. 

Knowing the potassium and argon content given weight mineral 
and the decay constants the age the mineral can calculated. 


The Problem Argon Loss 

The argon measurement will correct the following three conditions are 
satisfied. 

(1) None the argon has escaped before the sample was heated the 
container. 

(2) All the argon was extracted from the mineral. 

(3) None the argon was lost during the purification process. 


The last these conditions easily checked carrying out yield runs. 
sample was placed the container, run the normal way, and the volume 
argon measured. This measured volume was then taken back the U-trap 
and frozen down again. second identical sample was placed the container 
and also run the normal way. After purification the measured volume was 
twice the volume the single sample. This showed conclusively that none 
the sample was lost during the purification process. 

more difficult prove that all the argon being extracted from the 
mineral, but the experimental results indicate that this the case. The 
consistency the results obtained with the sodium hydroxide flux one such 
indication. This supported results obtained from the same samples run 
with solid sodium flux. The ratios obtained for the latter samples 
are also given, Table IV. this case the amount argon extracted increases 
the sample becomes more finely crushed, and the volume extracted from the 
finest sample closely approaches the mean result obtained with the sodium 
hydroxide flux. This indicates that even with solid sodium flux complete 
extraction being approached under these conditions. Potassium feldspars 
and lepidolites, different minerals with completely different lattice structures, 
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TABLE 
EFFECT VARYING SAMPLE MESH SIZE 


Solid sodium flux Sodium hydroxide flux 
size grains runs runs 
0.140 
60-100 0.021 0.070 0.137 
0.017 0.093 0.137 
200-250 0.010 0.108 
0.119 
250 0.002 0.134 0.137 


Mean standard deviation 0.138 .002 


both give ages agreeing with established uranium ages. This suggests that the 
lattice structures have both cases been completely broken down the 
sodium hydroxide flux. 

The good agreement between branching ratio values calculated from one 
our measurements and measurement Wasserburg and Hayden (48) 
constitutes check the whole extraction and purification method (40). 
The two laboratories use entirely different purification and measuring methods, 
and hence the agreement can most easily explained assuming that both 
laboratories are extracting and measuring all the argon contained the 
samples. Wasserburg and Hayden use sodium hydroxide flux and argon-38 
spike isotopic dilution method first used Suess, Hayden, and Inghram 
(45). used McLeod gauge measurement after complete purification 
the argon the manner already described. 

The last important indication that all the argon extracted the agreement 
obtained between the branching ratio measured counting experiments 
(Table II) and the results obtain using minerals known age (40). 
The mean result obtained from counting experiments 0.09 .038 and the 
result obtained from two minerals known age was 0.089. This agreement 
much better than could expected view the large scatter the 
results from the counting experiments; may purely fortuitous. 

has been suggested that part the argon may driven into the walls 
the calcium furnace glow discharge caused the induction heater. 
have operated three different calcium furnaces with two different induction 
heaters and have never observed any glow discharge. This possibility argon 
loss was checked the yield runs just described and also sparking the 
argon means Tesla coil. neither case was argon loss observed. 

Loss argon from the sample before its actual extraction may have occurred 
during the geologically long period since the formation the mineral. During 
that time argon could possibly have been lost diffusion into the surrounding 
rock. difficult test this experimentally but the consistency our 
experimental results for lepidolites and feldspars greatly differing ages 
indicates that any possible diffusion loss must small. Loss argon may also 
have occurred during the much shorter period from the crushing the sample 
and the placing the sample container the evacuation the system. 
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show that argon was lost owing the initial pumping different samples 
the same mineral were pumped for min., hr., and hr. The results 
agreed within the experimental errors. 

The possibility argon loss due the crushing was tested grinding 
samples different grain sizes. The results these tests are shown Table 
IV. can seen that the amount argon extracted was constant within 
the experimental error when the sodium hydroxide flux was used. This indicates 
that argon was lost from the sample because the crushing. 


RESULTS 


Six ages obtained the method described above are 
shown Table Every sample was run least twice and the standard 
deviation the individual argon measurements given. Where other ages 
existed for the same localities, these are shown the table 

The Lee Lake perthite was obtained from pegmatite few miles north 
the northeast shore Lac Ronge. uraninite was obtained from the same 
pegmatite. This uraninite has been dated Toronto the lead-ratio method 
and age 1740 million years obtained. The age 
for this sample 1750 million years. This sample one two used 
recently determine the branching ratio potassium-40 (40). can seen 
.that the decay constants which are present using give 
argon age excellent agreement with the lead-ratio value. 

The Silver Leaf Mine small pegmatite mine the Winnipeg River 
few miles upstream from the power dam Pointe Bois. has various 
times been called the Bear Claim and the Bob Claim and shown under 
these names provincial claim maps. situated approximately one-half 
mile from Huron Claim, which has been dated many times the various 
methods. these determinations only those based 
the radiogenic ratio have given consistent value. uraninite 
and monazite from the Huron Claim were analyzed isotopically Nier 
(33). From his figures calculated ages 2475 and 2590 million years 
respectively for these minerals. second uraninite obtained independently 
the writers has been analyzed this laboratory and lead-ratio age 
2590 million years obtained. The age 2550 million years 
for the Silver Leaf Mine lepidolite appears excellent agreement with 
the lead ages for Huron Claim. 

Holmes (23) has recently published the following age determinations made 
Rhodesian monazites: 


207 206 207 208 
206 
Ebonite Claims, Bikita, 
Southern Rhodesia 2680 2675 2680 2645 
Jack Tin Claims, north Salisbury, 
Southern Rhodesia 2650 2260 2470 1940 


Irumi Hills, Northern Rhodesia 2620 2040 2330 1340 


i 
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generally the case for old samples, the ratio gives the only 
consistent ages. The excellent agreement for ages the Bikita sample suggest 
that this very precise determination. Although this locality some 
miles from Popes Claim, there some reason for believing that the areas may 
parts single orogenic belt (23). single sample microlite has been 
obtained from Popes Claim and lead-ratio age determined for Toronto. 
The age obtained was 2510 100 million years. The age 
for Popes Claim lepidolite was found 2650 million years. Again the 
agreement between the two methods geological age determination seems 
very satisfactory. 

the writers’ knowledge uranium minerals are available from the 
Yellowknife area for age determinations. However number ages have been 
calculated from galenas found this area the method suggested Russell, 
Farquhar, Cumming, and Wilson (36). The mean eight these age deter- 
minations 2160 million years with standard deviation 150 million years. 
The age obtained the method for the Moose Dike peg- 
matite 2160 million years. The method appears best 
advantage such areas this where uranium minerals are exceedingly 
rare. 

Dill Township, Ontario, the Grenville Geological province. The age 
the Grenville orogeny has been assumed between 800 and 1100 million 

(7). The age 900 million years therefore does not 
seem unreasonable. 


CONCLUSIONS 


The measurements this paper have been presented illustrate the 
enormous potentialities the method dating and show 
that present techniques are adequate provide dates with fair precision. 
not suggested that the method yet established, nor that the ages 
obtained will always correct. However, wherever comparisons can made 
between ages and ages calculated other methods, the 
agreement striking. 

The rate beta decay for potassium-40 was determined averaging the 
results number counting experiments. The branching ratio used was 
previously chosen give agreement between lead-ratio ages and 
argon ages for two samples analyzed independently two laboratories. 
The fact that these same decay constants give good agreement for samples 
other ages indicates that they may well chosen. More precise counting 
results would provide much more stringent test the 
method. 

The assessment ages complicated the fact that 
various age determination methods often give discordant results for the same 
location. The writers know method that gives invariably good results, 
nor for that matter, any method that gives invariably bad results. However 
good estimate the correct age can often made from number deter- 
minations carried out the same area. Such comparisons have shown that 
for ages over one thousand million years the lead-ratio method gives much 
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more consistent results than the other uranium thorium methods (7). 
For this reason prefer use lead-ratio ages date our old samples wherever 
such ages exist. The “‘galena (36) less well established, but has 
given apparently good ages for old lead minerals, and have used these 
ages for the Yellowknife area where other ages have been determined. 

studies are being continued this laboratory, and 
several others. seems certain that many new determinations will reported 
very soon and that further tests the method will pos- 
sible. also likely that the method will soon applied potassium- 
bearing rocks well minerals. Such tests may have 
important bearing the future course geochronology. 
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H-PLANE BIFURCATION RECTANGULAR 


ABSTRACT 


Using method based the calculus residues, rigorous solution has been 
obtained for the problem the bifurcation rectangular waveguide. Expres- 
sions are given for the amplitudes all the reflected and transmitted modes 
the guide. comparison made with results obtained the transform method 
Wiener and Hopf. 


INTRODUCTION 


the solution electromagnetic boundary value problems all but the 
most trivial type, infinite sets equations infinity unknowns may 
arise. The solution such sets often very difficult, and usually approximate 
methods have resorted to. However, recent years technique for the 
exact solution certain class these equations has been developed 7). 
This method, which based upon the calculus residues the theory the 
complex integral calculus, seems offer powerful tool for the solution 
many such problems. 

The problem considered here the bifurcation rectangular guide 
semi-infinite septum parallel the electric vector the incident 
modes. This problem has been solved rigorously the method 
(4). The method used here makes use the calculus residues and some- 
what simpler and more straightforward. 


STATEMENT THE PROBLEM 


The configuration studied shown Fig. rectangular guide 
width fitted with perfectly conducting, infinitely thin septum starting 


Fic. H-plane bifurcation waveguide. 


from some reference plane and extending infinity towards the right. 
For the purpose analysis the space inside the guide will divided into 
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three subregions (see Fig. 1): subregion represents the space inside the 
main guide for subregions and represent the spaces the two 

assumed that with the electric vector parallel the 
edge the septum incident the plane from each the three 
subregions. (This assumption somewhat more general than usually 
necessary practice.) system reflected modes will set the three 
regions, with amplitudes and phases adjusted satisfy all boundary condi- 
tions. Since the incident fields are independent the y-coordinate and the 
entire structure cylindrical with respect the y-axis, the total fields the 
three regions will also independent The only field components present 
will and H,. The entire field can derived from two-dimensional 
scalar wave function 

The problem, therefore, find solutions the two-dimensional wave 
equation 

subject the following conditions: 

(a) and are finite all three subregions, except the edge the 

(b) and are continuous all three subregions and across the interface 
atz 

(c) represents outgoing waves apart from the incident 
components. 

(e) goes zero the edge the septum approached the 
distance from the edge). This implies that the magnetic field) goes 
infinity the edge (2, 5). 

ANALYSIS 

The following wave functions satisfy conditions (a), (c), and (d) regions 


n=1 
C sin b) > C, in b) 
n=1 
The first term each expression represents the fields, and the 
and are the propagation constants the three regions. They are given 
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where the phase constant for free space.* most practical cases, 
where all propagation constants will real except for and 
perhaps Hence, one mode only can propagate region and one 
mode may propagate region but not both. One the amplitudes 

Continuity the tangential field components across the interface 
(condition (b)) requires that 


p=1 


p=1 


n=1 
p=1 
where 
[5] 
forn 


One may consider the expressions the right Fourier sine series 
representing periodic functions period and respectively, which are 
ively. The coefficients these sine series are determined the usual manner: 


noted that the series are not absolutely convergent, since 
the coefficients A,, B,, C,, will seen later, behave for large 
The validity the term-by-term integration involved obtaining the 
Fourier coefficients may therefore questioned. 

possible, however, show that Equations [6] are correct. Green’s 
second identity, 


*Time dependence assumed. 
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while the contributions from and vanish when the radius the 


Jl 


Fic. Contour integration. 


semicircle and tend zero. The integrations along DE, involve only 
the absolutely converging series Equations [2], and yield Equations [6a] 
and [6c] the limit similar argument leads Equations and 
[6d]. 

deriving Equations [6], has been assumed that not rational. 
When rational, some terms the series will become indeterminate. 
easily seen, however, that Equations [6] and all the following equations 
(including the final results) are still correct, these indeterminate terms are 
replaced their limits a/b approaches rational number. 

When and are eliminated from Equations [6], two sets equations 
are obtained for the coefficients viz.: 


2 


Solutions must found for the which satisfy Equations [7] simultaneously 
and which lead fields that satisfy the edge condition (e) the first section. 
Similarly, when and are eliminated from Equations [6] one obtains 


[8] pr Bp n=1 On+ Bp 


. 
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SOLUTION THE EQUATIONS 
Consider function f(w) complex variable which analytic except for 


Cauchy’s theorem, 


f(w) can selected that the integrals the left tend zero 


then Equations [8] become formally identical with [7] and one obtains 
[11] 


The vanishing the contour integrals [9] equivalent the 
are excluded. 

find explicit expressions for and one starts with the contour 
integrals 


These contour integrals vanish Because [10] and [11], Equations 
[12] are then equivalent Equations [8] one sets 


seen that once the function f(w) known, the coefficients A,, B,, 
and can found explicitly from [11] and [13]. obtain the right poles 
and zeros, must have the form 


where g(w) entire function, the infinite product 


[13] 


[14] 
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and II(w, and are similar infinite products with and replacing 
and and replacing The factors and the exponential factors 
have been introduced make the products absolutely and uniformly conver- 
gent. 

manner similar that used Whitehead (7). found that 


the neighborhood the positive real axis. ensure that f(w) 
both halves the w-plane, excluding the neighborhood the 
poles f(w), g(w) must the form 

where P(w) polynomial of, the most, third degree. will shown 


below that P(w) must, fact, second-degree polynomial that the only 
remaining boundary condition, the edge condition (e) the first section, 


satisfied. 


The formal solution the problem thus given 
2 
eoW +e,wteo II(w, B) II (w, 7) 

and the Equations [11] and [13]. The coefficients eo, and can found 
from the three linear equations obtained substituting [19] [10]. 

The asymptotic behavior |w| given 


The second expression gives the correct poles and zeros the positive real 
axis. 
INVESTIGATION THE EDGE CONDITION 

establish that the above solution unique, necessary show that 
the electric field goes zero that the magnetic field goes 
infinity the edge the septum approached. The second condition 
easier verify. equivalent the requirement that the current density 
edge. 

From [2] and [13] the current density the side the septum 


q 
q 
{ 
{ 
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series differs only bounded function from the series 


[22] 
n=1 


zero considering the inequality 


[23] 


The integral the right the one the left 

The proof for the current density the other side the septum 
analogous. Similarly, one can show that has the correct singu- 
larity the edge when the approach made from region this case, the 
residues for large are found from the second expression [20]: 


From [2] and [11] then seen that differs only bounded function 
and from the series 


1 


DISCUSSION SPECIFIC CASES 


practice the guide dimensions will such that purely imaginary 
region but all other a’s are real. One may then assume without loss 
generality that all the for and all the y’s are real. Only two cases 
where the guide wavelength 

will assumed that unit amplitude incident from 
region but waves are incident from the right 
From [10] and [19] one obtains 


and from [11] the reflection coefficient referred the plane 
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where the bracketed superscript the product sign means that the factor 

Case real. From [27] and [15] seen that |R| i.e., there 
complete reflection, could expected from physical reasoning since 
waves are propagated regions and The junction introduces phase 
shift equal 


R=A,= 


’ ’ 0) +5,(22 ’ a ’ 0 ’ 
where 
n=N (n ) n 


This equivalent removing the septum and placing short circuit 
the guide agreement with the result the Waveguide 
Handbook (Ref. 174). 

Case II: purely imaginary. this case one obtains 


Similar expressions are obtained when wave incident from the right 
circuit and agrees with that given the Waveguide Handbook (4, 305). 


[30] 
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ABSTRACT 


have developed method increasing the concentration He* 
mixtures “superfluid relying upon the fact that liquid He‘ flows rap- 
idly through small leak while the solution remains behind. this way 
increase concentration from 95% can achieved single operation. 
processing about 200 standard cc. the initial gaseous mixture 
per hour. 


INTRODUCTION 


Enrichment the content mixtures has been achieved 
the conventional methods thermal diffusion (2, 13, 15, 16) and fractional 
distillation (7). However, the failure liquid become superfluid like 
liquid He‘ suggested entirely new methods. The method “superfluid 
(5, was based the observation that liquid flows rapidly through 
narrow channel, while liquid does not. The method (7, 11, 
14, 17, 19) took advantage the fact that the moves with the normal 
component the liquid away from source heat. 

most these investigations, however, the initial mixture was either 
atmospheric helium containing part well helium containing 
about part 107 and the final concentration never exceeded 
one case (14) the thermal diffusion method was first used obtain 
concentration 0.01% and the heat flush method was then applied achieve 
final concentration 4%. The work described the present paper started 
with mixtures produced by-product the operation the Chalk River 
nuclear reactor, which contained between and Our task was 
the concentration least 25% order that the mixtures might 
useful for nuclear bombardment experiments, and also hoped that 
might achieve still higher concentrations because the interest attached 
the properties these mixtures the liquid state. chose the method 
superfluid filtration because Abraham, Weinstock, and Osborne (1), the 
course their work the mixtures, had obtained 
appreciable enrichment this method the concentration range that 
interested us. practice the method proved simple and effective, and 
were able develop the point obtaining concentration 95% after 
single filtration. 

short preliminary account the earlier stages this work has already 
been published (3). detailed account has been given the Ph.D. thesis 
Lovejoy (12). 

Throughout this paper, percentage concentration defined 100 
where and are the numbers gram moles He* and 
present. 
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APPARATUS AND METHOD 


The apparatus shown diagrammatically Fig. The superleak was 
the type originally used Giauque, Stout, and Barieau (9), and was 
made sealing cm. length 0.18 mm. diameter platinum wire into 
pyrex glass. Upon cooling low temperatures the platinum shrank away 


Vacuum 


Separation Receivin 
Cryostat Cryosta 


perleak 


Fic. apparatus. 


from the glass, leaving gap which was found have mean width about 
cm. observing the rate flow helium gas through liquid 
helium temperatures. the final apparatus such superleaks were used 
parallel. The platinum wire was wound lengthwise short length pyrex 
tube which was then placed inside another pyrex tube, filled with pyrex 
rod, and the whole assembly melted down and drawn out. 

measured quantity low concentration mixture from the container 
was liquefied top the superleak, which was immersed liquid helium 
bath 1.3°K. Flow through the leak started immediately, but ceased soon 
the pressure the exit side became comparable with the saturated vapor 
pressure. The explanation this may that the liquid mixture the exit 
side had lower concentration than the entrance side, and the resulting 
osmotic pressure (4, 20) was large enough prevent flow. When the exit 


Vacuum 
Helium 
Line 


J 
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side was pumped and the pressure there kept below the saturated vapor 
pressure prevent the formation bulk liquid, the flow started again. 
the amount the effluent was sufficient justify its recovery, 
ordinary pumping system could not used, and was not possible 
design Toepler pump with high enough speed. The tube from the exit 
side the superleak was therefore taken out the separation cryostat 
into the receiving cryostat which was maintained the slightly lower 
temperature 1.1°K. The effluent therefore distilled over from 
and collected liquid the tube Cryostat was equivalent pump 
with high speed limited only the diameter the connecting tube. the 
end the run the impoverished mixture was boiled off into the container 

Although the liquid flowing through the superleak contained little 
its concentration was much less than the initial mixture, and the con- 
centration the liquid left behind the entrance the superleak 
increased steadily. The vapor pressure mixture increases with 
increasing concentration and has been measured concentrations 
13% (18). was therefore possible follow the course the enrichment 
observing the vapor pressure the manometer Curve Fig. 
shows that the concentration increased slowly first, then more rapidly, 


Vapor pressure mm. 


Time minutes 


Fic. the vapor pressure the liquid the entrance the superleak. Curve 
the amount liquid which has passed through the superleak, deduced from the height 
the liquid column the tube The initial charge was 250 standard cc. mixture. 
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and finally levelled off the flow came end. Curve shows the amount 
liquid which passed through the superleak function time, and was 
deduced from the level the liquid the tube The manometer enabled 
check kept the concentration the effluent. 

When the filtration had come end and the vapor pressure registered 
had reached its maximum value, the enriched mixture was removed 
the Toepler pump into the container Its concentration was measured 
with mass spectrometer. 


PERFORMANCE 


The maximum vapor pressure attained the end filtration always 
corresponded liquid concentration 2%, but the measured con- 
centration the enriched mixture was always considerably greater than this. 
This easily understood once realize that the vapor equilibrium with 
24% liquid mixture has concentration about 95% (18), and that the 
dead space above the superleak liquid helium temperatures was the 
order 0.75 and therefore held considerable mass vapor. follows 
that there optimum value for the quantity gas initially introduced 
one wishes realize the maximum final concentration. smaller quantity 
gas introduced, all the liquid flows through the superleak before liquid 
concentration 24% reached. larger quantity than the optimum 


used, then liquid remains the end the filtration and lowers the mean 


concentration below that the vapor. For our apparatus, determined 
experimentally that this optimum quantity was 200 standard cc. 
mixture 100 standard cc. mixture. 

However, even when this optimum initial charge was used, the resulting 
mean concentration was still not 95% but about 45%. This anomaly was 
traced the fact that the gas the dead space room temperature had 
concentration the range having experienced appreciable enrich- 
ment. not yet clear whether this was because low concentration mixture 
room temperature can remain equilibrium with high concentration 
mixture liquid helium temperatures, because the rate diffusion 
from the low the high temperature was very slow. The difficulty was sur- 
mounted skimming off the gas room temperature and then removing 
the cold gas separately. was this way that final concentration 95% 
was achieved. For the procedure which gave this high concentration, the rate 
processing was 200 standard cc. the initial mixture per hour, yielding 
about standard cc. 95% mixture per hour. Only 50% the the 
initial charge was contained the 95% mixture, 25% having been lost through 
the superleak and 25% removed with the gas room temperature which had 
been skimmed off. 

attempt was made increase the concentration above 95% this 
method. probably more feasible bridge the gap between 95% and 100% 
fractional distillation. have been able raise the concentration from 
80% 92% merely liquefying the 80% mixture glass capillary. 
1.1°K. and pumping off the enriched vapor with Toepler pump. Starting 
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with high concentration mixture, Fairbank, Ard, Dehmelt, Gordy, and 
Williams (8) have used fractional distillation produce concentrations better 
than 99%. 

The concentration the effluent collected the receiving cryostat depended 
upon the initial concentration and also the initial charge complicated 
way (12). This impoverished mixture could also processed and the final 
concentration was found almost independent the initial concentration. 
The rate processing was, fact, slightly higher, but, course, the yield 
cubic centimeters high concentration mixture per hour was much less and 
the procedure became unduly laborious when the initial concentration was 
less than 0.3%. 

The fact that was able flow through the superleak should not neces- 
sarily taken evidence that can participate superfluid flow. The 
possibility ordinary viscous flow the cannot excluded. The situation 
may have been complicated thermal gradients and concentration gradients, 
and the effective pressure forcing the through the leak may have been 
quite large. Moreover, there some evidence (10) that the effective viscosity 
the normal component very small channels this width, perhaps 
because mean free path effects, and similar considerations may apply 
solution. 
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THERMAL EXPANSION LITHIUM, 77° 300° 


PEARSON 


ABSTRACT 


Measurements the lattice parameters lithium temperatures below 
300° have been made and the results are compared with previous work and 
considered terms the Griineisen relationship. Methods low temperature 
X-ray photography are also discussed. 


INTRODUCTION 


Accurate measurement expansion coefficients very low temperatures, 
say K., matter some difficulty and present not many data are 
available, although there considerable interest studying the validity 
the Griineisen relationship low temperatures (see, for instance, Bijl and 
Pullan (5)), and obtaining the expansion the inert gases. 
X-ray cameras have been built this laboratory measure accurately 
lattice parameters low temperatures and particular obtain the ex- 
pansion coefficients the alkali metals part general program carried 
out towards understanding the electronic nature these ‘simple’ metals. 

The present paper describes lattice parameter measurements lithium 
the temperature range and discusses methods used low tempera- 
ture X-ray photography this laboratory. 


LOW TEMPERATURE X-RAY METHODS 


The camera used this investigation between 80° and 300° consists 
collimator system and base plate which film cassette precisely located, 
adjustable specimen mount, and synchronous motor for rotating the 
specimen. The X-ray film held cm. Unicam high temperature camera 
film cassette first employed Basinski and Christian (4). The camera can 
used either with the Post, Schwartz, and Fankuchen (17) Hume-Rothery 
and Strawbridge (8) method cooling, while temperatures 77° 90° can 
obtained running liquid nitrogen oxygen over the specimen after the 
method Smith and Lonsdale (12). The specimen its adjustable holder 
may either hung supported rampant, being rotated bearing accord- 
ingly either above below the X-ray beam. 

method suggested Dr. Dugdale has also been used. The X-ray 
camera placed large cylindrical can having cellophane windows for 
entry and exit the X-ray beam, Fig. Although specimen icing effectively 
prevented this method, suffers from the difficulty that the camera 
cassette cools somewhat below the temperature which its dimensions have 
been measured. Trouble may consequently result through differential con- 
traction the aluminum alloy film cassette and the steel knife edges. This 
may change the standard angle between the high angle knife edges and 
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~THERMOCOUPLE LEADS 


GLASS TUBE 


SPECIMEN 
RUBBER CORK 
CM. CASSETTE 
ELLOPHANE COVER 
CELLOPHANE COVER 
BEAM IN 
SPECIMEN 


SPECIMEN CENTERING 
STAGE 


MOTOR LEADS 


MOTOR 
Fic. arrangement for taking low temperature X-ray powder photographs. 


also the standard circumferential distance between the low angle knife edges. 
Variation not normally allowed for taking account film shrinkage 
but Lipson and Wilson (10) have shown that the errors involved are negligibly 
small. Changes might indeed expected lead inaccuracies lattice 
parameter determination. However, calculations embracing both changes 
and relative camera diameter show that the error the extrapolated 
lattice parameter obtained from reflections high Bragg angles negligibly 
small for any probable degree cooling the film holder—especially 
relation the uncertainties involved precise determination specimen 
temperature. 

order measure the temperature which photographs are made, 
touch the specimen just below the level the X-ray beam. This thermo- 
couple may either calibrated fixed points, situ measuring the 
lattice parameter substance whose contraction already well known 
the relevant temperature range. The second method preferable for number 
reasons, among which may mention: 

(1) Calibration the thermocouple takes place under the exact temperature 
gradients obtaining exposure any particular temperature. 


— 
He 
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(2) Any number exposures and subsequent calibration may made 
relate the characteristics particular thermocouple the standard e.m.f. 
tables used. 

(3) Any differences between thermocouple temperature 
specimen temperature the particular temperature the exposure are taken 
into account. 

values Nix and MacNair (16) obtained for 99.99% aluminum 
intervals approximately have been used calculate the change 
lattice parameter aluminum with temperature. (The few earlier values 
obtained Ebert (7) and used Hume-Rothery and Strawbridge (8) are 
approximate agreement with those Nix and MacNair.) Calculations have 
powder 99.99% aluminum kindly supplied The Aluminum Company 
Canada Limited. this way calibration curves have been established while 
using the various methods cooling. guard against errors arising from the 
temperature calibrations, have made measurements lithium using two 
different methods cooling: (1) the arrangement shown Fig. and (2) 
modified form the Hume-Rothery and Strawbridge type cooler. The 
agreement both sets measurements, after using the temperature cor- 
rections obtained from the appropriate calibration curves, certainly well 
within the absolute experimental errors temperature measurement, which 
are probably less than +4° the range 

this work lattice parameters are given for ready comparison with 
earlier work. have used radiation for all photographs (Acux,; 
1.537395 kX). X-ray films have been measured with Cambridge universal 
measuring machine using low magnification, and lattice parameters have been 
calculated using the Nelson and Riley (15) extrapolation function, applied 
resolved doublets occurring high Bragg angles. The lithium (240) 
Lattice parameter measurements are accurate +0.0002 for the best 
films and perhaps +0.0004 for the remainder. Combined with absolute 
temperature uncertainty say +4° mentioned above (1° 0.00016 
expansion for lithium room temperature) encouraging find the 
results two series experiments different methods self-consistent 
+0.0002 over the whole range. 

The lithium used was special sodium-free metal prepared New Metals 
and Alloys Limited and used Kelly and MacDonald (9) for resistance 
measurements and MacDonald and Pearson (13) for measurement 
thermoelectric force. The residual resistance the metal, 
1.5 and spectroscopic analysis the actual X-ray specimens used 
shows the presence calcium and very slight traces sodium the only 
metallic impurities. order prepare X-ray specimens, the metal was 
extruded into vaseline prevent oxidation and lengths were sealed thin- 
walled pyrex capillaries which were evacuated. The metal was annealed either 
for long time room temperature. Several specimens were 
prepared. One these used the range 140-300° showed considerable 


surface discoloration. 
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RESULTS AND DISCUSSION 


Fig. are shown the lattice parameters two specimens measured 
different techniques the range 125-300° K., together with that another 
specimen measured 77.5° lattice parameter 3.5016 kX, corrected 
20° C., obtained for lithium. This lower than the value 3.5023 


3-5000 


100 200 300 
OK. 


Fic. Low temperature variation lithium lattice parameters 
measured values ‘a’ using two different methods cooling. 
calculated variation ‘a’ with (See Text.) 


found Lonsdale and Hume-Rothery (11), but agrees with the estimate 
3.5015 3.5019 for pure sodium-free lithium made Aruja and Per- 
litz (1). The difference the present result and that Lonsdale and Hume- 
Rothery, which lies the limit the probable experimental errors, appears 
due difference purity the lithium samples.* 

Lonsdale and Hume-Rothery measured the lattice parameter lithium 
90° K., and the observed contraction this temperature, 0.0261 
kX, agrees with the present findings, 0.0258 kX, measure- 
ment one part ~12,000, which satisfactory. Only approximate agree- 
ment found with the results macroscopic measurements made Simon 

*Resistivity measurements and spectroscopic analysis show that the lithium used this investi- 


gation virtually sodium free. Less than 0.1% sodium solid solution the lithium used 
Lonsdale and Hume-Rothery would account for the parameter difference. 
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and Bergmann (18). also noted that the expansion coefficients copper, 
iron, and nickel, determined Simon and Bergmann the same method 
used for lithium and reported the same paper, are not agreement with 
subsequent work, e.g., see comments Nix and MacNair (16). 

The present measurements may further assessed terms the Griineisen 
relationship The temperature variation the linear thermal 
expansion coefficient shown Table obtained for assumed value 


TABLE 
Calculated Calculated 
Temperature, linear expansion Temperature, linear expansion 

300 46.8 180 36.8 
290 46.4 170 35.6 
280 45.7 160 34.4 
270 45.0 150 33.0 
260 44.3 140 31.4 
250 43.5 130 29.6 
240 42.7 120 27.7 
230 41.8 110 25.4 
220 40.9 100 23.0 
210 40.0 20.1 
200 39.0 17.3 
190 38.0 


0.92 taking the values for lithium calculated from the measurements 
Simon and Swain (19) Kelly and MacDonald (9), values the 
volume occupied gm. molecule, from the present work, and values 
the compressibility, calculated from Bridgman’s (6) measurements 75° and 
30° C., assuming (in the absence other data) linear temperature variation. 
The line plotted Fig. calculated from these expansion coefficients, 
normalized value 3.5019 295° K., and shows satisfactory agree- 
ment with the measured lattice parameters over the range 77° 300° 
The value 0.92 considerably lower than 1.17 given Mott and 
Jones (14). Further discussion not possible until the exact temperature 
dependence has been established for the temperature range investigated. 

Simon and Bergmann (18) obtained very good agreement calculated 
and measured expansion coefficients using the Griineisen relationship and 
values corresponding 510° (Simon and Swain’s (19) specific heat 
measurements lead almost constant value 400° 15° for 
lithium temperatures from 100° 300° (Kelly and MacDonald (9).) 
Kelly and MacDonald have pointed out that lithium anomalous that 
above 100° (resistance) smaller than (specific heat). The inclusion 
Schottky transition leading higher intrinsic value, suggested 
Simon and Bergmann result their expansion coefficient measurements, 
would rather accentuate the anomaly since presumably Schottky 
transition would not affect the electrical resistance. 

order examine further the reported Martensitic transition (Barrett and 
Trautz (3)), b.c.c. Li, two photographs were taken tempera- 


| 
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tures approaching that liquid helium running stream liquid helium 
from syphon directly over the rotating specimen using the apparatus shown 
Fig. The helium flow was adjusted maintain thermocouple touching 
the specimen just below the X-ray beam axis constant temperature which 
was the lowest that could obtained matter how fast the helium flow was 
made. The temperature the specimen was believed not greatly above 

The main fact emerge from this experiment was that the b.c.c. form was 
still present large extent these temperatures. The lattice parameter 
the b.c.c. form 3.472; .002 (Fig. 2). Only limited 
information could gained about the h.c.p. phase because the probable 
coincidence reflections (102), (104), (203, 120), (114, 122), etc., with those 
the b.c.c. form, and (100), (101,002), (200, 103) with those and possibly 
LiOH either these were present impurities. However, its presence could 
certainly inferred from the (110) and (201, 112) reflections medium and 
weak intensity respectively and 36.6° which were free from the 
possibility overlaps. Similar relative displacements the interplanar 
spacings the h.c.p. phase those recorded Barrett and Trautz (3) are 
found and these have been shown Barrett (2) due the presence 
stacking faults. 

could not definitely established whether there was any the h.c.p. 
phase present specimens photographed 77.5° 
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RADIATION FROM SLOT 
CYLINDRICALLY TIPPED 


JAMEs WAIT AND SIDNEY KAHANA 


ABSTRACT 


solution outlined for the electromagnetic field distribution magnetic 
current and charge volume bounded infinite wedge with cylindrical 
tip both infinite conductivity. The result applied calculate the fields 
thin half-wave transverse slot the cylindrical surface. Radiation patterns 
the equatorial plane are calculated for wedge angles 90°, 180°, and 270° for 
various radii the cylindrical surface. The results have application the design 
slotted cylinder antennas for directional coverage. 


INTRODUCTION 


Extensive use now being made the radiating properties slots cut 
the surface metal cylinder. Slotted cylinder antennas are particularly 
useful very high and ultra-high frequencies obtain broad coverage (5, 
6). For the particular case axial slot, Sinclair (9) has shown that the 
radiation pattern the horizontal essentially circular the diameter 
the cylinder less than about quarter wavelength. The pattern 
becomes more directive the diameter the cylinder increased. means 
thereby exists which the horizontal pattern can controlled some 
extent. Sinclair also shows that variety patterns can obtained 
feeding power different amounts pair axial slots the periphery 
the cylinder. 

alternative method, suggested Alford (1), control the radiation 
pattern slotted cylinders employ metal wings the slots. The main 
effect this shaped reflector intercept the current which normally 
flows around the cylinder and force into the wings. The horizontal pattern 
then rendered more directive with the radiation being confined mainly 
the region subtended the 

the purpose this paper derive expressions for the radiation 
arbitrary slot with specified tangential voltage, the cylindrical surface 
which coaxial with infinite corner reflector wedge. This idealized 
theoretical model the winged slots mentioned above. The solution this 
problem may also interest connection with slots radiating from 
airfoil and tail surfaces aircraft. The final results are quite similar the 
expressions derived Papas (5) and Silver and Saunders (8) which are 
restricted slots isolated cylinders. 


THE FORMAL SOLUTION 


The general problem calculate the fields (see Fig. due toa 
magnetic current element with respect conducting wedge angle 
with cylindrical tip radius Choosing cylindrical polar coordinate 
system, the surfaces the wedge are defined and the tip 
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Fic. Wedge with cylindrical tip. 


and are (p, and ¢’, 2’) respectively. The solution the problem 
most readily carried out Green’s function method. Expressions for the 
field vectors are then obtained terms the currents and charges the 
region and the fields the surfaces the region. this particular problem 
more convenient deal with magnetic currents and fictitious magnetic 
charges rather than their more conventional electric counterparts. 

Following Stratton (10) readily shown that the electric and magnetic 
field vectors E(r) and H(r), within region bounded surface 
(with outward normal are given 


and 


where J*(r’) and p*(r’) are the magnetic current density vector and magnetic 
charge the point P’. The operator involves differentiation with respect 
the primed coordinates. The Green’s function G(r, r’) satisfies the inhomo- 
geneous Helmholtz equation 


where the three-dimensional Dirac function (3). has been assumed 
the above that the time factor exp(iwt) and the medium contained 
the volume homogeneous with electrical constants and 

simplify the form equations [1] and [2] convenient impose 
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certain boundary conditions r’)/dn over the surface 
For example the axial component the magnetic field can written 


subject the conditions p=a, and 
therefore seen vanish, which result the assumed perfect conductivity 
the wedge and the vanishing contribution from the closing spherical 
surface infinity. 

With similar reasoning follows that 


The essential task this stage find the forms the Green’s function 
which satisfy the above-mentioned boundary conditions. can seen 
that G(r, r’) can written the following general form 


1 +00 


Fourier coefficient G,, seen satisfy the following inhomogeneous equation 


which for points Bessel’s equation order Solutions equa- 
tion [7] which are continuous give rise outgoing waves infinity, 


(2) 


order and the Hankel function the second kind order The 
primes the Bessel function indicate differentiation with respect the 
argument. The corresponding form G,, satisfy the boundary condition 
Bessel functions are removed. The factor can found integrating with 
respect both sides equation [7] over small interval which includes 
This process yields C,, The formal solution the problem 
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then obtained when the required expressions for the Fourier coefficients are 
substituted into equation [6]. 

particular application this result calculate the fields distribu- 
tion magnetic current the surface the cylinder. Schelkunoff’s 
equivalence principle (7) such distribution equivalent specified tan- 
gential electric field the surface the cylinder. The process obtain this 
special result straightforward remembered that equation 
continuity relates the magnetic charge and the magnetic current vector. 
also convenient make use the Wronksian relation (2) 


when the transition made from the Green’s function. Omitting 
further details the derivation follows from equations [4] and [5] that 


+co 

m=0 —oo 

and 
+o 

m=0 


The coefficients a,, and b,, which are functions are related the prescribed 
electric field components and the portion the 
forms these coefficients are given 


The remaining field components can derived the same process they 
can obtained directly from equations [11] and [12] noted that the 
field can always represented (11) basic set (transverse magnetic) 
and (transverse electric) waves. They are given 


m=0 
+o) 

m=0 
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The integrations with respect can now carried out under the usual 
far field approximation where terms order 1/R? and higher are neglected. 
The results this saddle-point integration (12), expressed terms spherical 


polar coordinates (R, are 


Inserting the values the coefficients, and the radiation fields are 
given 


and 
z oa(2’) 
; iv cos 0 ikcos@2’ 


apparent that the component the radiation field depends only 
the axial electric field E,(a, the slot. also seen that this axial 
slot voltage also produces transverse component the radiation field. 
the other hand the transverse electric field the slot produces 
only component the radiation field. 


THIN TRANSVERSE SLOT 


The above general result can now specialized thin slot 
cylinder (see Fig. 2). can assumed that the slot sufficiently thin 
the voltage distribution varies sinusoidally along the slot, that 


o— gw 


for slot centered where the voltage the center the 
slot. further consequence the assumption very thin slot can 
seen that the contribution from the transverse voltage which proportional 
vanishingly small. The radiation fields the thin half-wave 
transverse slot centered the bisector the wedge and located 
the cylindrical surface are then given 
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where 

noted that the equatorial plane (i.e. 90°) the component 
vanishes leaving only the component the field. The radiation pattern 
proportional this plane called the principle pattern the trans- 


Fic. Pattern for transverse slot isolated cylinder. 


P 
z 
| 
120° 
180 A 
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verse slot. The following cases have been calculated: and and wedge 
angles and These results along with those for isolated 
cylinder are illustrated Figs. 


Fic. Pattern for 270°. 


The shape the patterns rather interesting. When the cylinder isolated 
with wedge, the field behind the cylinder seen considerably reduced 
relative the forward field. This effect much more pronounced than for 
the case axial slots where the patterns are almost circular the equatorial 
plane (8, 9). would expected the diffraction effects become more pro- 
nounced for the larger circular cylinders. 

The patterns for the slot the cylinder with the wedge the reflector 
present are quite different. these cases noted that the field always 
zero directions along the surface the wedge, which compatible with 


270° 
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the boundary condition that the tangential electric field zero these 
surfaces. important feature these patterns for the wedge that for 
larger values ka, multiple lobes can occur, which are usually undesirable. 

The rigorous extension the above theory wedge reflector finite 
extent would very difficult indeed, not impossible, present-day 
techniques. believed, however, that the assumption infinite reflector 
does not limit the practical application the results much first sight. 
The effect diffraction around the extremities the reflector somewhat 
similar the case vertical electric antenna parallel the apex 
corner reflector (13) where has been observed experimentally (4) that the 
magnitude the spurious field behind the reflector within few per cent 
the central maximum the reflector made square sheets dimen- 
sions greater than two wavelengths. fact, consideration knife edge 
diffraction would predict that the field behind the reflector the equatorial 
plane would the order times the central field, where the 
length side the sheet. 
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NOTE 


EXPERIMENTAL STUDY BAND INTENSITIES THE FIRST POSITIVE 
SYSTEM 


III. QUANTITATIVE TREATMENT EYE ESTIMATES 


NICHOLLS 


INTRODUCTION 


the purpose this note examine what quantitative information 
may obtained molecular constants and source conditions from eye 
estimates band intensities. Naturally, for the complete interpretation 
molecular spectrum terms these quantities, accurate intensity measure- 
ments the bands are primary requisite. Such measurements are however 
seldom available, and yet quite often eye estimates band intensities have 
been reported together with vibrational analyses band systems (5). 

spite the recognized deficiencies subjective estimates intensities, 
seems necessary examine whether any quantitative information may 
obtained from these estimates when other data intensities are available. 
This situation often occurs the preliminary interpretation low dispersion 
spectra weak sources. 

example, some eye estimates intensities the first positive system 
have been considered and have been shown yield substantially the 
same results, far the behavior electronic transition moment with 
internuclear separation concerned, were obtained parts (7) and 
(8). these papers, accurately measured band intensities were used together 
with Factors’ manner suggested Fraser (2) find 
the dependence electronic transition moment upon internuclear separation. 


GENERAL CONCEPTS 
The integrated intensity band may written 


where constant depending upon the optical system and units employed, 


and the other symbols have been previously defined and II. 

the photographic plate used working the roughly linear region 
its characteristic, the photographic density produced the band may 
written as** 


*This terminology for the overlap integral square has been suggested Bates (1). The vibrational 


the band sufficiently wide that the band head more simple measure than the whole 
band, may defined the band head density, and then will replaced equation 
the intensity included the band head. For bands one system, may expected 
roughly constant from band band and dependent rotational temperature (4). 
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where the contrast factor the emulsion response, and the inertia 
the emulsion wavelength Further, the action the eye estimating 
blackening compared that the photocell microdensitometer, 
application the Weber—Fechner law sensory responses stimuli suggests 
that the rough relationship existing between estimated blackening and photo- 
graphic density 

const+CD,, 


where constant which determines the scale (usually 10) adopted 
estimation the blackening. 
Thus from equations [1], [2], and [3] 


where Cy, modified contrast factor. For the purposes the following 
discussion will assumed that adequate data (3), (5), and 
(2) are available. 

METHOD 


Equation [4] indicates that the following simple procedure may reasonably 
used determine, albeit somewhat crudely, the transition and 
the source from eye estimates band intensities. 


Criteria Quantitative Significance 

The subjective nature the eye estimate, the possible variation plate 
response with wavelength, and also, particularly small dispersions, the 
possible unsuspected presence underlying band structure other systems 
make very necessary the establishment some criteria quantitative signifi- 
cance applied the data before use. 

Equation [4] suggests that these criteria may applied connection with 
plate response were roughly constant over the wavelength range bands 
progression, each plot would expected linear, slope and 
relatively spaced ordinate appropriate the influence the log term. 
Scatter points about the straight line for one progression would due 
the influence the log term. Slow variation plate sensitivity will 
tend disturb the rough equality slope from line line while really 
unreliable estimates intensity and also the presence unsuspected under- 
lying band structure will cause marked deviation individual points from the 
lines their progessions. This also useful way uncovering, from 
intensity point view, unsuspected underlying bands. 

Data for aberrant points that lie far from their respective lines should 
course rejected untrustworthy. regions where plate response rela- 
tively constant across the wavelength range the system, average slope 
can then found from the slopes the plotted lines. 


Determination the Electronic Transition Moment 
For the remaining data, plots should made 
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for each the progressions. From equation [6] will seen that this 
equivalent plotting 


for each progression over the range appropriate it. The plots for each 
progression will vary relative ordinate owing the log term, and 
order show how log varies with over the range encountered 
the whole consistent displacement ordinate must applied 
each the plots allow for the effect the log term, and make 
possible for smooth curve drawn through the points displaced. this 
way the behavior log with may found, and only cases where the 
resulting plot line parallel the axis one led believe that the electronic 
transition moment probably independent internuclear separation. 


Determination Vibrational Population 

cases when examination the degree vibrational excitation the 
source which consistent with the eye estimates necessary, the following 
method may used. From the scaled plot log versus described 
above, value log for each band may read from the smooth 
curve. According equation [4] plots 


should, for progressions, linear, slope unity, and displaced ordinate 
from the line Thus the vibrational population may 
found from these displacements and plotted versus v’. Such 
plots will show, their tendency to, departure from linearity, whether 
not the vibrational excitation roughly form and whether 
certain vibrational levels have been favored population depopulation. 


EXAMPLE 


While the above treatment has been applied with success number 
band systems whose eye estimates are available, perhaps the best test its 
valid use seen the comparison between the dependence upon 
for the first positive system found from accurately measured intensities 
described and and that found from eye estimates intensity. 

The eye estimates have been taken from Pearse and Gaydon’s well-known 
compilation data molecular spectra (5). The estimates there recorded for 
this system are from two sources: (a) Poetker’s (6) band head intensities 
greater than about 7000 measured thermoelectrically, and (b) Pearse and 
Gaydon’s estimates blackening Profs. Fowler and Strutt’s 
plates less than about 7000 A). When the test significance described above 
was applied, showed that these two sets numbers could not considered 
the same scale equally representative. Accordingly, the following 
discussion confined solely the eye estimates. 

For these data, application the test significance led the rejection 
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eye estimates the (4,0), (5,0), (6, 1), (7, 2), (8, 3), (10, 7), (11, 8), 
(12, 6), (12,9) bands quantitatively unreliable. For the 
remaining bands the plots log +log versus 
were made and the rescaling procedure was employed. was found that 
representative straight line consistent with 

[5] R.(r) const exp(—2.7r) 


could drawn through the rescaled data. will recalled that the compara- 
tive result obtained was 
const exp(—3.02 r). 


may also mentioned briefly that for the second positive system 
and for the first negative system analysis recently recorded integrated 
intensities (9) has led substantially the same R,(r) curves for these systems 
were obtained from analysis the eye estimate data (5). 

This, together with the substantial agreement between equations [5] and 
[6], indicates that the method described this note for treating eye estimates 
band intensities can used some cases obtain roughly the same 
trends with may obtained when more accurate intensities are 
available. 

Details the application this procedure other band systems and 
the rough evaluation vibrational distributions will published shortly. 

This work was performed connection with Contract 
held with the Air Force Cambridge Research Centre. 
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